On the Detection of a Scalar Stochastic Background of Gravitational Waves 



Danilo Babusci 
INFN- Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 

Luca Baiotti 
SISSA and INFN sez. di Trieste, 
Via Beirut 4, 34013 Trieste, Italy 

Francesco Fucito 
INFN sez. di Roma 2, 
Via della Ricerca Scientifica, 00133 Roma, Italy 

Alessandro Nagar 
Dipartimento di Fisica, Universita di Parma 
and INFN, Gruppo Collegato di Parma, 43100 Parma, Italy 



o 

(N 
C 

o 

(N 



O 

in 
o 



In the near future we will witness the coming to a full operational regime of laser interferometers 
and resonant mass detectors of spherical shape. In this work we study the sensitivity of pairs of 
such gravitational wave detectors to a scalar stochastic background of gravitational waves. Our 
computations are carried out both for minimal and non minimal coupling of the scalar fields. 



I. INTRODUCTION 



(N ■ 

In few years, the research on the detection of gravitational waves (GWs) will hopefully greatly progress. This hope 
is based on the coming into operation of a new generation of experimental devices that, if they can be operated at the 
planned sensitivity, should probe deeply into the region in which we believe GWs can be observed. According to the 
experimental technique employed, these detectors can be divided into two categories: interferometric detectors and 
resonant mass detectors. To make our point more concrete let us concentrate on Michelson interferometers Jl],|2| and 
resonant mass detectors of spherical shape 0. The main advantage of interferometers is their sensitivity in a wide 
frequency band. On the other hand spherical shaped resonant mass detectors at resonance have the same sensitivity 
""^ ' regardless of the direction of the impinging GW. 

q-( In the following we will concentrate on a very specific issue, that is on the possibility of detecting scalar GWs. Our 
interest in this subject stems from the observation that Einstein's gravity is definitively not the only mathematically 
consistent theory of gravity and in fact the presence of scalar fields coupled to gravity is required by a vast array 
of theories that model various phenomena as the inflationary universe or attempt to incorporate gravity with the 
.J^ | quantum world. For a review on this subject see Ref. Q. For more recent proposals that also require a modification 
of Einstein's gravity see Refs. (5|-0|. Are all the above described detectors fit to measure scalar GWs? While the 
answer is obvious for resonant mass detector of spherical shape, the situation for interferometers must be analyzed 
with care. Let us use for a moment the "standard" description, that is well suited for our kind of argument, of an 
impinging GW (for the moment we neglect its spin content and direction) stretching the lengths L\ and L2 of the 
two arms of the interferometer. The conventional Michelson interferometer is configured for maximizing its sensitivity 
in the detection of the differential mode signal A_ = SL\ — SLi- Even if the information regarding L\ and L2 
separately is available, the sensitivity of these measurements is orders of magnitude worse than that of A_, and, thus, 
a single interferometer of this type is not able to disentangle the common mode signal A + = SLi + 8L2 (transverse 
monopole mode) from A_ (usual spin 2 mode). A way out could be the construction of an array of these detectors 
or the adoption of a different optical configuration (Fox-Smith) for the interferometer Even if interesting from a 
theoretical point of view, these alternatives do not seem practical, given the cost and the difficulty in operating such 
complex apparatus. A viable alternative to these proposals could be, from our point of view, that of a coincidence 
analysis on the data of an interferometer and a resonant mass detector of spherical shape || . 

In this work we study the sensitivity of combined pairs of resonant mass detectors and interferometers to a scalar 
stochastic background of gravitational waves (SBGW). If such a background has a flat spectrum (which is the standard 
assumption) even the narrow frequency band available to a resonant mass detector won't have much influence on our 
conclusions. Our computations generalize the results of Ref. || in which the sensitivity patterns to scalar radiation 
were considered. To be as general as possible, the impinging radiation is computed in the general setting given by 
scalar tensor theories . Our main result is the computation of the sensitivity to scalar GWs of correlated pairs of 
(solid mass or hollow) resonant mass detectors of spherical shape or pairs of interferometer-resonant mass detector. 
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Finally we consider the effects on such detectors of massless non minimally coupled scalar fields, generalizing the 
results of Refs. [ p~T|Jl^1 . While this paper was being written, a similar analysis employing two LIGO interferometers 
for massive and nonrelativistic scalar particles appeared p3]. 



II. SCALAR TENSOR THEORY 
A. Fundamental equations 

Let us consider a very general tensor multi scalar theory of gravity, where the gravitational interaction is mediated 
by n long range scalar fields (p a in addition to the usual tensor field present in Einstein's theory. The action in the 
Einstein frame is 

\-i / J4 



S=(16nG)- 1 J d i x^{R-2g^ lab {^)d^ a d vV b ) + S m [* m , A 2 {tp a ) 9liv ] . (2.1) 

We use units in which the speed of light is c = 1 and the signature is — +++. Greek indices A, /i, v, . . . = 0, 1, 2, 3 denote 
spacetime indices; Latin indices from the second part of the alphabet i,j,k,l... = 1, 2, 3 denote spatial indices; Latin 
indices from the first part of the alphabet a, 6, c, . . . = 1, . . . , n label the n scalar fields. Our curvature conventions 
follow those of Ref. R = g^ u R pl , is the curvature scalar of the Einstein metric <7 M „ and g = det(g (Llw ). The action 
contains a dimensionful constant G, which will be denoted as the bare gravitational constant (related to G Newton's 
constant as measured by Cavendish experiments) and a a model type metric 7 a b(^), not necessarily positive definite, 
in the n dimensional space of the scalar fields. S m denotes the matter action, which is a functional of some matter 
variables and of the Jordan-Fierz metric g^ v = A 2 (ip)g filJ . The scalar fields can be non minimally coupled to 
matter. This means that they can appear as coupling "constants" between the matter fields ty m and gravity g^ u . For 
instance, low energy string type theories naturally introduce in the action terms with couplings of the kind 

Sm = "f / tfvyfl'pF&F&P'F*, (2-2) 

where = d^A^ - d v A^ + f ABC A^A^ is the Yang-Mills field strength and the scalar field <p is the dilaton. 

By varying the action S with respect to the Einstein metric g^ u and the scalar fields <p a , one obtains the following 
field equations 

~ \ Rg^u = 2 7a6 (v>) (d^ a d v ip b - \g^g pa d P v a d a ^ + 8ttGT^, (2.3) 
gTVf7 v <p a +g^jU<p)9^ b 9 u ip c = -47rG[a a (p)T + a a ] : (2.4) 
where 7^, are the Christoffel symbols of the metric 7 a &(y>). The functions a a ((p) = d a In A(ip) represent the field depen- 

— 1/2 

dent couplings between scalar fields and matter within the metric sector of the theory. T^ v = 2 (—5) SSm/Sg^ is 
the stress energy tensor, T its trace and a a = (—g) 5S m /5(p a is the density of scalar charge. In the Jordan-Fierz 
frame we would have 

T„ v = A 2 { v )f^ Vl o a = A 4 (<^)a a , (2.5) 

as can easily be found from their definition fIo|| . Actually, since \J—g = A 4 (ip)^/— g, and 

S v S m = J ' ^x^a a 5 V a = J d i x^/^a a 6ip a , (2.6) 

one recovers immediately ( [2.5[ ). 

In the literature, scalar tensor theories with a a = (i.e. metric theories) have been studied by many authors, from 
the pioneering work of Jordan, Fierz, Brans, Dicke, Wagoner jlBj to the recent studies of Damour and Esposito-Farese 
fiof . This interest arises from the fact that they do not violate the Weak Equivalence Principle and so imply geodesic 
dynamics for neutral weakly self gravitating bodies. However this is not the most general framework, in particular it 
is not the case of the interesting scalar fields foreseen by string theory. For a recent analysis see Ref. [jllj . 

Let us compute the expression of the relative acceleration between two weakly self gravitating bodies in the general 
n scalar theory; this formula will be the starting point to write the response of a GW detector to a scalar tensor wave. 
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When a a ^ 0, the stress energy conservation law in Einstein units is |T( 

V„T^ = a a V^ip a T - cr a W° , 

or, in the Jordan-Fierz frame 

V%„ + a a X7^ a = 0. 



(2.7) 



(2.8) 



This equation implies a non geodesic motion of test mass bodies. This result corresponds, for a single scalar field and 
a particular choice of the coupling function A(ip), to the lowest order gravidilaton effective action of string theory 
JyJ . However, if a a — we have 'V u T flv = and so geodesic motion of test mass bodies is recovered. 

In Ref . , starting from the single field string like case of ( |2.8| ) the equation of motion of test mass bodies has 
been derived. Following the same line of reasoning, we generalize that result to our case. Let us recall the point like 
limit of the generally covariant energy momentum tensor for a particle of mass m and world line x^ (r) fl4| 



P V~9 



(2.9) 



where p M = mdx^ /dr. We can rewrite the scalar charge density a a for a test body, in terms of dimensionless scalar 
functions q a , which express the relative strengths of non universal scalar to tensor forces 



v a {x') = -q a T(x') 



Qa ■ 



*< 3 >( 



5 (r)). 



(2.10) 



As we con sider long range fields, q a -C 1 to avoid conflicts with the present test of the Weak Equivalence Principle. 
From (2.8) we get the geodesic equation in scalar tensor theory with non minimal couplings [pT| 



T^ v x a x v + q a d^ a = 0, 



(2.11) 



where x^ = da;' 1 /dr. Now we can compute the modifications to the relative acceleration between two test mass bodies 
moving along two worldlines induced by the q a 's. 

Let us take two weakly self gravitating bodies moving along two infinitesimally close worldlines x^(t) and 
x'^(t) — x^(t) +S^(t), where <5 M is the separation vector between the two curves. If we suppose that the bodies 
have different scalar couplings qtP and , their relative acceleration is Q 



R 



IOJO 



(2.12) 



where Si = dSi/dt. Notice that in (2.12) there is a term proportional to q^P — q y a' ■ This term will be important when 
the test mass bodies are of different nature (e.g. one is a baryon and the other one a lepton) but it is irrelevant inside 
a GW detector. Therefore, the equation needed to analyze the response of GW detectors to scalar tensor waves is 



-(2) 



Rio jo + Qadidjip 11 



(2.13) 



B. Gravitational waves 

Let us recall some results concerning scalar tensor GWs |l(| . In the weak field limit of the theory 

/(x)=^ + f(x), (2.14) 

where \h^\ <C 1, |£ a | <C 1, r/^ is the flat Minkowski metric, and the background values of the scalar fields. We 
now choose a gauge in which the metric perturbation has zero time-time and time-space components while the purely 
spatial components, for a plane wave propagating along the direction characterized by the unit vector Cl, assume the 
form 

fc«(x) = h A (x)e? j (n) + 2ale(x)e s ij {£iy, A = +,x; a = l,...,n. (2.15) 
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e + , e x are the spin 2 polarization tensors describing the ordinary GW in the transverse traceless gauge, e s is the 
spin polarization tensor of the scalar waves, a° a = a a {tf>o), and we choose units such that A(ip$) = 1. By indicating 
with rh and n a pair of orthonormal vectors lying in the plane perpendicular to O, these polarization tensors can be 
written as follows (see Appendix |C 2) 



(2.16) 



and 



e^(f2)e s ' y '(r!) = 26 



BB' 



B = +, x, s 



We consider now the small relative oscillations of two weakly self gravitating bodies induced by t his wave. By 
indicating with Li the rest separation of the bodies, we can put Si = Li + Q <C 1). Expanding ( 2.13| ) to first order 
in (i, we find 



d 2 hj_ 
dt 2 



u 



(2.17) 



Since we are considering plane wave solutions, the spatial derivatives appearing in the last eq uatio n can be replaced 
by the time derivatives, namely didjt; a = (li£lj£ a = (Sij — ^-(f)))^, and taking into account ( 2.15|) , one finds 



ld^ 

2 dt 2 



h A (x) ef J (tl) + 2t 



q a )C(x)etAn) + 2q a C(x)Si 



and then the infinitesimal displacement induced by the GW is 



h A {x) e£(n) + 2 (a° a - q a ) FixWyifi) + 2q a t a {x)5 K 



(2.18) 



(2.19) 



This formula needs a few comments. The scalar fields considered in our theory are massless, therefore the scalar 
GW can carry energy and momentum thr ough just one degree of freedom, the transverse polarization tensor efj (fi) 
(see Wagoner in Ref. |lj|). Therefore in (2. IE) only the transverse part strains the matter and the Sij is effectively 
unimportant when studying the response the antennas to GWs. By introducing the effective gravitational wave sensed 
by the test mass bodies 



■eff 



h A (x) 4(0) + 2 (a° a - q a ) fWef/O) 



we rewrite (2.19) as follows 



1 - 
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(2.20) 



(2.21) 



However, if the scalar fields were slightly massive, ther e wou ld be also a longitudinal polarization along the propagation 
direction of the GW and we couldn't drop the Sij in ( 2.1£ ). This scenario has been analyzed in |) ll-lj|, but in the 



following we won't consider it and just restrict our study to h 



eff 



III. INTERFEROMETERS AND RESONANT MASS SPHERICAL DETECTORS 
A. Response function of an interferometer to scalar GWs 



Let us consider a M ichelson type laser interferometer with two orthogonal arms of the same nominal length L\ = 
L2 = L. From (2.21), the signal at the output port of the interferometer (the strain of the differential mode) is 
proportional to the difference in the two path lengths, Ci — C2, induced by the wave and can be written in the form 

ra 



4 



h eff =htJ f V», (3.1) 

where I? is a traceless and symmetric tensor describing the geometry of the interferometer Q In the interferometer 
frame, namely the one where the corner station stands at the origin of coordinates and the x and y axes lie along the 
arms, this tensor writes 



V= - -1 . (3.2) 




The effective strain sensed by the interferometer is then split in a spin 2 and a spin part, proportional to the 
difference oP a — q a ; we can take explicitly into account the dependence of the strain from the angles (Q, </>) defining the 
direction VL of the incoming wave by introducing the angular pattern functions of the interferometer 



and writing the strain as 



F A (n) = e$(n)x>« F s (n) = ,;,.:<>rP'<. (3.3) 

h^f = h A (x) F A (Cl) + 2(a° a - q a ) F s ({l) . (3.4) 

B. Cross section for resonant spheres in scalar tensor theory 

We discuss now the cross section of a resonant sphere in the general scalar tensor theory. For spin 2 waves this 
result was obtained in Ref. J0| (see also Ref. [fl9| ). In recent years, this kind of detector (both solid and hollow) has 
been extensively studied as device able to analyze the spin content of GWs (see Refs. pp|-p2|). The calculation of its 
scattering cross section in the framework of the Brans-Dicke theory was carried out in Refs. p3|-p5||. The extension 
of the results of Ref. ]23j] to the general scalar tensor theory with minimal coupling is straightforward [^6| and in 
the following subsections 1 and 2, we will just sketch the steps and quote the results. Furthermore, in subsection 3 
we will repeat the calculations for the even more general case of q a ^ 0. For the sake of generality the direction of 
propagation of the wave and the antenna frame (defined by a triad of orthonormal vectors (x,y, zj) will be taken to 
be distin ct. The direction f2 = (9, (f>) of the incoming wave is identified by the relative orientation of the triad defined 



in ( 2.1 6| ) with respect to ( 



in = cos <p x + sin <p y, 

n = — sin</>cos# x + cos(f>cos9y + sin#£, (3-5) 
ft = sin 4> sin Ox — cos <f> sin 9 y + cos 9 z. 



1. Tensor GWs 



Consider a superposition of spin 2 plane GWs with wave vector and amplitudes Iia impinging on a spherical 
GWs detector 



= e^e*'*" + c.c. = h A e A u e lk " x " + c.c, A = +,x. (3.6) 



Note that hereafter e A are the polarization tensors written in the detector frame (x, y, z) (see Appendix C 2 for their 
explicit expressions). As usual we will use the so called quadrupole approximation, i.e. we suppose that the detector 
is much smaller than the wavelength of the impinging GW, so that only the first terms (quadrupole, for the tensor 



1 (3.1) is valid in the regime in which the wavelength of the impinging scalar GW is much bigger than the length of the arms 
of the interferometer. Given the resonant frequencies of our resonant mass detectors, this will be always the case in the present 
paper. For a more detailed discussion of this point see Jl7| ]. 
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component; monopole and quadrupole for the scalar one) have to be considered. Analogously to Ref. |l4| we find the 
expressions for the spin 2 scattering and total energy cross sections 



128ttG 2 



tot 



/ 



(3.7) 
(3.8) 



where a 2 , 



and fij = Tij(£l, /) is the (traceless) Fourier transform of the variation induced in the stress energy 



tensor of the sphere by the impinging GW|. Furthermore we will study resonant scattering, i.e. we will assume that 
the detector scatters only the impinging GWs with frequency / around the resonant frequency of one of its natural 
vibrational modes. This leads to a relation between a s ^ at and crjj°' and to another between fy and the sphere mode 



tensors 



= 7(/)e 



(3.9) 
(3.10) 



where r\ is the fraction of the total oscillation energy dissipated through emission of GWs (it can be calculated as a 
function of the detector internal parameters) and j(f) gives the frequency de pen dence of fy. The function j(f) is 
chosen so that the response of the antenna is resonant in frequency (see Refs. po| , p3f ). If the oscillation of the mode 
with angular momentum I — 2 has proper frequency]^] f n 2 and a bandwidth A f n2 , we find 

7(/)« . , j .. f=. (3.11) 

/ ~ Jn2 + «A/ n2 /2 



Substituting (|3.7| ) and ( p.ll[ ) respectively into (|c 
becomes 



) and ( 3.10 ) and combining the results, the total energy cross section 



o"h(/; n,l = 2) 



1 GMv 2 F„ 



2tt 



(/ - fra)' 



+ A L 



,/4' 



(3.12) 



where M is the sphere mass, v the velocity of sound in the material the sphere is made of, F n a constant depending 
only on the quadrupolar mode under scrutiny and on the sphere parameters (radius, density, material) [E3j, and 
G = (1 + ag)G is the effective Newton's constant measured in Cavendish like experiments. The results of Ref. |2^] 
in Brans-Dicke theory are recovered by setting «q = (2u>bd + 3) _1 where u>bd is the Brans-Dicke parameter. If 

^ (e) (defined m 



in performing this calculation we expand the tensors in the numerator of cr^°* in the detector basis T>\ 



Appendix CI), we find that the total cross section can be written as the sum of five terms, the total cross sections 
for any single vibrational mode of the sphere 



{ (/) 



e = 0, lc, Is, 2c, 2s, 



(3.13) 



and 



1 GMv 2 F 71 



2tt 



A 



L,2 



F^h A 



(/ - / n2 f + A 2 /re2 /4 Y. A >\ h A- 



(3.14) 



Z5^e^'(f2), explicitly written in Table | 



The angular dependence is enclosed in the pattern functions, F A 

For later purposes we will need also the integrated cross section. By integrating ( 3.12| ) we get 

E h (n; I = 2) = -J-j GMv 2 F n 
l + ai 



(3.15) 



2 In principle the expression of fij could contain also a term proportional to 2?^°' cx 5%j p^| , accounting for the trace of 
the polarization tensor (monopole excitation). But, since the trace en vanishes, in this tensorial part such a term gives no 
contribution. 

3 The index n = 1, . . . ,oo labels different solutions for the sphere eigenmodes with fixed angular momentum I El. 
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2. Minimally coupled scalar GWs 



The scattering and total cross section for the minimally coupled scalar part of a GW are 

to&al [ I m ? jj^j 

s " 5 ^(f,/)e(f,/) ' (3 - 16) 

where t; a (x,f) is the Fourier transform of the impinging scalar GW. We now have to decompose the scalar GW 
polarization tensor e*j in a quadrupole and a monopole part, as they excite different modes in the detector. The way 

to do this is by expressing e*j in the basis defined by the five real symmetric tensors , plus £>fy , proportional 
to the ide ntity tensor, because these tensors are directly related to the angular momentum of the excitation (see 
Appendix ^T] and Ref. 

Assuming again resonant scattering and noting that the resonance frequencies of the quadrupole and the monopole 
modes need not be equal, we have now two expressions for the variation of the stress energy tensor of the detector: the 
first, labelled fij (/; I = 0) and valid for the sphere monopole mode, is proportional to T>f^ and has resonance frequency 
/ = fnol the second, (/; I — 2), is proportional^ to J2e F^V^ and has resonance frequency / = f n % ^ f n o 

(/;' = o) = /?'(/Kr(/)€ 0) , (3.18) 



T, 



'J I 



%(/; I = 2) = «Kr(/)£ F i e):D i e) , (3.19) 

e 

where 0'(f) ^ (3"(f) are the analogous of the function j(f) in ( |3 . 1 OP and = T>lfe^(Cl). We deduce then the 
total cross section of the monopole mode 

^ = 0) = "\ dMv2Hn (f . /4 , (3-20) 

l + c*o 77 (/ - /nor + A/ n0 /4 

where i?„ is a constant depending on the monopolar mode under exam |23| ] and A/ re0 is the resonance bandwidth. 

ie calculation gives 

a s (/;rM = 2) = £V| e )(/,r2) (3.21) 



For the quadrupole modes , the same calculation gives 



where 



if - f n2 r + A}j4 

The pattern functions Fg are listed in Table S: since an explicit computation yields 



^ (3-23) 



the global response to scalar waves of the quadrupole modes is isotropic too, and total cross section ( 3.21 ) reads 

., (/;M = 2 ) = ^M»> F ft. (3.24) 

As the quadrupole modes are sensitive to scalar and to tensor waves, the angular dependence of each cross section 
could make it possible, in principle, to guess the polarization. For instance, considering the m = mode, Fs (9,<f>) 
gets a maximum for 9 = <p = 0, while f} 0) (0, 0) = F { ° } (0, 0) = 0. 



4 The I = 2 part of e? i expanded in the basis is 2j] e F^vf) . 
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TABLE I. Angular dependence of the sphere pattern functions for the three independent polarizations of a scala r te nsor 
GW. Notice that the pattern functions of the e — 2c mode coincide with the ones of the interferometer introduced in (3.3) 



Mode (6) F[ e \e^) Ff\8A) F. W (M) 

2s -cos#cos20 -i (l + cos 2 9) sin 20 - 1 sin 2<f> sin 2 9 

2c — cos 8 sin 2<f> \ (3 + cos 29) cos 2cj> i cos 2(f> sin 2 9 

Is — sin S sin i sin 2$ cos — i cos sin 2$ 

lc — sin 9 cos <f> — | sin 26* sin <fr ^ sin sin 29 

^ sin 2 9 ^(3cos 2 0-l) 



S s ( n; I = 0) = GMv 2 H n 



The integration of (3.20) and (3.24) gives, respectively 
and 

GMv 2 F n 



E s (n; Z = 2) = 



l + a 2 3 



(3.25) 



(3.26) 



3. Non minimally coupled scalar GWs: q a ^ 

In this section we present the full generalization of the result presented before to the case in which q a is small but 
not exactly null. We will follow step by step the procedure outlined in Ref . |^3) . Further details for the general multi 
scalar metric theory can be found in Ref. p6| . 

a. The energy momentum conservation law 

Fi rst l et us consider the energy momentum tensor T^ u of the resonant sphere and write the linearized conservation 
law in momentum space. Denoting by ?^ v (x), f (x) the linear part of T^ v ,T we get 

d^ u (x) + a a (x)d u e(x)=0, (3.27) 

which in momentum space reads 

Vr^(fc) + ?„(*) * [k"C(k)} = 0. (3.28) 

The reality of f ^jx) and £, a (x) implies f^*(k) = f^(-k) and £, a *(k) = £, a (-k), with k = (k,uj) and u = k = \k\. 
The asterisk in (3.28) stands for the four dimensional convolution product. Now we proceed to express fbo(fc) in terms 

Offy(fc). 

For a particle of mass m, (2.10) defines the relation between the scalar charge densities and the components of the 
energy momentum tensors. Integration over all particles of the resonant sphere gives 

5-„(a;) = -q a f(x), (3.29) 

and therefore the four equations ( |3.28] ) in momentum space read (with r M „(fc) = f^v) 

fc r 00 + k t f° l - q a f * [k°C{k)] = 0, (3.30) 

k f 0i + kjf ij - q a f * [fc^ a (fc)] = 0. (3.31) 

The wave trav els al ong the directi on Cl , and so fej = k^VLi because the scalar fields are massless. Subtracting the 
contraction of ( 3.31 ) with Cli from ( 3.3C ) gives then 

f 00 = %f!W, (3.32) 

a relation which holds in minimally coupled scalar tensor theories too p(| . 

We now compute again, in the case q a ^ 0, the quantities entering the cross sections: the incoming energy flux, 
the power emitted by the detector in GWs and the interference power (see Ref. |lj] ) . The calculation strictly follows 
that of Ref. for Brans-Dicke theory which has been generalized in Ref. |2(| to multi scalar metric theory. These 
latter results are recovered in the limit q a — > 0. 
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b. The incoming energy flux 

Let us start with the incoming energy flux which is independent from the direction of the incoming GW thanks to 
the symmetry of the detector. We will simplify here the calculations assuming the incoming direction to be coincident 
with the z axis of the detector frame. Later we will recover the general expression for an arbitrary direction. 

The incoming flux is computed gi ven t he energy momentum pseudotensor of the gravitational field. At second order 
in the linear expansion, defined in ( 2.14 ) 



t 



(2) = _±_ 

8ttG 



(daa b U a d^ b + -a° b ^(d p h vp + dji pp - dph^e + 



-(d a a b )or) l *»r 1 f " T Z a d p d„Z b + -agj^V" {d~h lp + d^h ep - d p h ie )d^ b + 



-[(d a a b )o + a a a b ~ 1 ab \d fi ed„e 



-a„a.i 



(d a a b ) - -^a" a a b - -%,, 



8nG 



(3.33) 



where 



R$ = \h ap ( dpdjiap - dp,d p h aiJ - dpdjiup + d a d p h pv 



+ | ( dphap + d a hpp - d p h ap ) ( d a h p + d v h ap - d p h% ) + 
1 

~4 



dp,K p + dvh^ - d p h MV 2d a h ap - d"h 



(3.34) 



is the Ricci tensor linearized to second order in the fields. We keep in mind that 

oo- 

£+ + 2a% a £ x 

£ x -£++2a°J a 

,0 0, 



(3.35) 



and denote by < ... > the integrat ion o ver a 3 dim ensional space region with linear dimensions much bigger than 
the GWs wavelength. Substituting ( |3.35| ) into ( 3.33| ) we obtain the total scalar tensor energy flux coming from the z 
direction 



(3.36) 



c. The scattering amplitude and the energy cross sections 

Let us consider a GW impinging onto our spherical resonant detector. At large distances, R — \x\, from the detector 



C(x,t) 



,,2-KifR 



R 



-2-wift 



(3.37) 



where A a (x,f) is the scattering amplitude relative t o th e ath scalar field. It obeys the usual reality condition 
A°(x, /) = A a *(x, — /). Using the scalar field equation (2.4), under the hypothesis that the quadrupole approximation 
holds, the scattering amplitude can be written in terms of f pl/ (x) as 

A a (x,uj)~G [ d 3 x' {a^-q a )f(x',uj)e- lls ' 
= Gl 



(3.38) 



where we have expressed fpo(fc ) in terms of the space like components of the Fourier transform of the energy momentum 
tensor by making use of ( 3.32| ). 

Let us turn then to the detailed calculation of the energy cross section, referring ourselves again to Ref . |2^] . From 
( |3.38 ) we find the scattering power to be 







P scat = ^ / dfiA (£,/)A-(x, /) = 



(3.39) 



where g 2 = g a g a - Furthermore, the interference between the incident plane wave and the scattered wave gives 



pint _ 2/ q. 



dtt£ Q (f,/)A Q *(x,/)<5(l-fc-x) 



47T/S {C&f) (aS-g")f«e« 



(3.40) 



The scattering and total cross sections are then 



_„,.„,. _ P scat _ 8-kG- ( o 



(a 2 , - 2g Q ag + g 2 ) 



{ir«r+!^f«} 



pint 2G Q {K-<?a)e*%e^} 



/ 



(3.41) 



(3.42) 



where we have put £ c = £ c (ir, /). Expanding now in the Vf>°\ X>« basis, we can decompose into an I = part 
and an ^ = 2 part 



W;J = 0) = C(/) K-ga)r^° 0) ' 

nj (/; I = 2) = C"(/) (a° - g.) C £ F i e) . 



/ — fnO, 
f = fn2, 



(3.43) 
(3.44) 



with C'(/) 7^ C"(/) defined as /?'(/) and /?"(/) in ( 3.18 ) and ( [3.19 ). Hence, the monopole and the quadrupole total 
cross sections become 



a s (f; n,l = 0)= — 9f(C ) ^ - 



2G 



ffa (/;n,/ = 2)=- r 9f(0 E( F « 



2 L(e;a2,g a ) 



where 



l (r ; a°, ? a ) = i a° Q r i 2 + 1 q a c i 2 - K % r e + & «g r* ^ 



(3.45) 
(3.46) 

(3.47) 



By using (3.41) and ( |3.43[ )-( p.46| ) with the analogous of ( |3.9| ) with a s replacing a^, and assuming once again resonant 
scattering, we get the final form for the monopole and quadrupole total cross sections 



a s (f;n,l = 0) 



a s (f;n,l = 2) 



A L/ 4 L(e;*Ua) 



irp (a 2 - 2g a ag + <f) (/ - /„ ) 2 + & ^ 



15 7/2 



^/ 2 (a 2 - 2g a ag + g 2 ) (/ - /„ 2 ) 2 + A 2 n2 /4 



E(^ (e) ) 



(3.48) 
(3.49) 



We still have to evaluate 770 and 772- This is done remembering their definition as the ratio between the power 
pscat — p(n;l) reem itted as gravitational waves by the vibrations of the sphere and the oscillatory energy E^lp 
dissipated by the sphere itself 



To 



p(n;0) 



p(n;2) 



2%A fn0 E' o 



(n;0) ' 



'72 



2tt A fn2 E' 



(n;2) ' 



(3.50) 



The oscillatory energy is that evaluated in Refs. [^0|j23j, since it doesn't depend on q a . The calculation of the reemitted 
power follows that of Ref. |2^,2(|. The only difference consists in replacing a° a with a° a — q a . Therefore, omitting the 
uninteresting details of the calculation, we get 
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770 = 4G {a - 2q a a + q }, 



% = 



A/„ 
2G Mv 2 f 2 2 F n 



{a 2 Q -2q a a% + q 2 }. 



15 A /nL 

Finally the cross sections assume the following simple forms 

GMv 2 H n A fn0 L(Z a ;a° a ,q a ) 



v s (f;ri,l = 0) = 



(/ - / n0 ) 2 + A}J4 



<r.(f;n,l = 2) 



GMv 2 F n 



A, 



In (f-f n2 ) 2 + A 2 j4 



2 L(e;ag,gq) 



(3.51) 
(3.52) 

(3.53) 
(3.54) 



These expressions can be made more manageable by expanding L g a ) in powers of g Q <C a" -C 1, an ordering 

relation which follows from the weak field limit of (2.7). 

First, an analogous calculation to that for A a (x, /) gives [ p6| 

G 



C(x,f) = -(a a -q a )n i e l J, 



(3.55) 



where now = fL(k) is the Fourier transform of the space components of the stress energy tensor of the source 
located at a great distance R from the antenna; f ^ is related to the Fourier transform of the variation of the quadrupole 
moment QL (/) of the source by 



Q«(/) = - 



2tt 2 / 2 v ' 

Therefore, taking into account ( 3.56| ), an explicit evaluation gives 



.an -a 4tt 4 / 4 G 2 



L(£°; tt 2>9«) 



i? 2 



|<&e?| 2 a$ 1-4 



q a ao 



,q a a^q 2 q 4 



and 



C« e = ±L ^-|0 / «e?l a ogh-2 



qa a$ 



Expanding this ratio in powers of g a yields 

L(?;a° a ,q a ) 



ag - 2 q a a a a + q 2 + 



(3.56) 



(3.57) 



(3.58) 



(3.59) 



We can finally compute the q a dependent terms in the cross sections. 
The monopole cross section at the lowest order in q a reads 



<r s (f;n,l = 0) 



1 GMv 2 H n 



A 



/no 



1 + a 2 



(/ - /„ ) 2 + Ajj4 



{a 2 -2q a a a } 



(3.60) 



where we have reintroduced the effective Newton's gravitational constant G. 
Analogously, the quadrupole cross section for any mode e writes, at first order 



1 GMv 2 F„ 



1 + al 



2tt 



{f-f n2 ) 2 + A 2 jA 



Summing over e, (3.23) gives 

= 2) 



1 GMv 2 F n 



1 + ag 



67T 



(/-/ n2 ) 2 + A 2 n2 /4 



{a§-2? a8} 



(3.61) 



(3.62) 
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IV. DETECTION OF A STOCHASTIC GW BACKGROUND 



Our aim is to generalize the standard analysis about the detectabili ty o f the spin 2 stochastic GW background 



[p8| [3 Of to the case of the general scalar tensor theory outlined in Sec. II B . Within this framework we introduce a 
density of scalar gravitational radiation p s in addition to the standard tensor one pu ■ If we assume, as in the tensor 
case, that the scalar background is isotropic, unpolarized, stationary and Gaussian, it is completely described in terms 
of the (dimensionless) spectrum 

n - = 1 dT7 ■ ^ 

Pc d In / 

where dp s is the energy density of the scalar gravitational radiation in the frequency range / to / + df and p c is the 
critical density required (today) to close the universe 

Pc = ■ (4.2) 

8irG V ' 

Ho is the present value of the Hubble constant. Notice that, although we study a scalar tensor theory we normalize 
the scalar gravitational spectrum to the value of p c recovered in general relativity. This choice has been taken to 
have a direct comparison between the tensor only and the scalar tensor framework. The present value of the Hubble 
expansion rate is usually written as Hq = ho x 100 km s^ 1 MPc^ 1 , where ho (= 0.6 4- 0.7) is a dimensionless factor 
that parametrizes the experimental uncertainty affecting the value of Ho- As a consequence of this definition the 
quantity /iqQ s (/) is independent of ho, and, thus, more suitable to characterize the stochastic GW background. 



A. The signal to noise ratio for scalar tensor GW stochastic background 

From the experimental side, the signal induced in the detector output by a stochastic GW background is indistin- 
guishable from the intrinsic noise of the detector itself. Unless the amplitude of the signal is very large, then, the 
subtraction of an a priori estimate of the detector noise cannot be confidently applied to the data. This implies that 
in order to detect a stochastic GW background, we should rather analyze the correlated fluctuations of the outputs 
of, at least, two detectors with no common sources of noise (a condition usually verified for widely separated detector 
sites). The cross correlation among detectors is advantageous also from the point of view of the minimum detectable 
signal. It can be shown |2^ , |3l|| that, under the same experimental conditions, the minimum detectable signal in the 
correlation of two detectors can be even three orders of magnitude smaller than the one detectable with a single 
detector. 

The problem of the optimal processing of the the detector outputs for the detection of the stochastic GW background 
(tensor and scalar) has been considered by various authors pp8|-|30|] , and extensively reviewed in Ref. fnl. This 
analysis can be generalized with minor modifications to the case of the general scalar tensor theory considered here. 

The signal present at the output of each detector can be written as (we consider the case of two detectors) 

s k (t) = n k (t) + h e k ff (t), (4.3) 

where we have indicated with h e J^ the gravitational strain due to the stochastic GW background and with n the 
intrinsic noise of the detector, while k — 1,2 labels the detector to which each quantity is referred. The noise 
is assumed to be stationary, Gaussian and statistically independent on the gravitational strain. Furthermore, the 
assumption that the noises in the two detectors are uncorrelated implies that the ensemble average of their Fourier 
components satisfies 

<nl(f)n l (f)>= S(f-f)6 kl ^Si k \\f\), (4.4) 

where S^\\f\) is the (one sided) noise power spectrum for the kth detector. Given an observation time T, the 
correlation "signal" is defined as follows 

,T/2 ,T/2 

S= dt At' Sl (t)s 2 (t')Q{t-t'), (4.5) 

J-T/2 J-T/2 
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where Q is a real filter functionH that, for any form of the signal, is chosen in order to maximize the signal to noise 
ratio (SNR) associated to S [p0[ 7 



The statistical treatment of the signal S defined in (4.5) starts with the plane wave expansion of the metric 
perturbations. The effective GW exciting the detector at position x writes 



df / dQ 

Is 2 



2irif(t-€l-x) 



(4.6) 



where, as a consequence of the reality of h^j (t, x), we have that h* A (f) = hA{—f) and £ a *(/) = £ a (— /); is the unit 
vector specifying the direction of the incoming GW. 

In terms of the detector tensor, the GW strain sensed by the detector k located at Xk is given by 



hl ff (t) EE h'"(t,S k ) = h^ f (t,X k )V», 



and (4.7), keeping into account (4.6) and the definitions (3.3), it can be rewritten as 



df / dQ 

Is 2 



h A (f,n)F k A (n) + 2(a° a - qi k) )C(f^)Fi(n) 



c ,2irif(t-Q,-x k ) 



(4.7) 



(4.8) 



In the following, we focus on the spin contribution to the strain, i.e. the part of ( |4.8| ) depending on £ a (/, VL) (the 
spin 2 contribution has been extensively treated in Ref. J28|-|3(|). Since the scalar background is assumed stationary, 
Gaussian, isotropic and unpolarized in the space of the scalar fields, it can be shown (see Appendix [A] for further 
details) that the correlation functions between the Fourier amplitude £ a (/, O) of the waves are 



<r(/^K b (/',tt) > 



(4.9) 



where f\(/) is the spectrum of a single scalar field. As a consequence of our assumptions on the scalar background, 
the whole spectrum is then fi s = nfl^. Following the same line of reasoning applied in the case of the spin 2 waves, 
under the further assumptions that the detector noises are much larger in amplitude than the gravitational strain and 
statistically independent on the strain itself, for the SNR we obtain 



SNR ? = 



a Q - \ q a +q a J a Q + q a q b 7 3^2 



1 + al 



8tt 3 



«I(/)r|(/) 



1 1/2 



-'/ / df 

f*sP(f)sP(f) 



(4.10) 



The function T^(/) is the generalization to scalar fields of the usual overlap reduction function introduced in p9| , ^0[ . 
This is a dimensionless function describing the reduction in sensitivity due to the different location and orientation 
of the two detectors, and it is given by 



s 2 



dQ Fl(CL)F£(tt)e 



2rrif d£l-s 



(4.11) 



where s is the unit vector along the direction connecting the two detectors and d is their distance. Notice that T^(f) 
coincides with the scalar overlap reduction function introduced in Ref. jt^ in the context of the single scalar metric 
theory of Brans-Dicke. The SNR obtained in Ref. is also recovered specializing (4.10) by setting q^ = and 
a 2 a = {2u BD + 3r 1 . 

In the following, we will consider two detectors with the same qa — qi^ — q a - If we keep only linear terms in q a , 



(4.10) becomes 



SNR; 



al - 2g a qg 3ff 2 
1 + a 2 8tt 3 



2T 



df 



^ 2 (/)r|(/) 



M 13 (/)si 2) (/) 



1/2 



(4.12) 



5 This function depends only on t — t' as consequence of the assumed stationarity of both the gravitational strain and the 
detector noise. 
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As pointed out before, the spin 2 contribution to the strain ([D|) is substantially the one already obtained in the 
literature. However the presence of the scalar fields slightly modifies the usual formula for the SNR/j of Ref. |3(i| l 
introducing an overall a§ dependent factor which is absent in General Relativity 



SNR ft 



3#o 2 



8tt 3 



IT 



f 6 sL 1} (f)si 2 \f) 



1/2 



Here &h{f) is the usual spin 2 spectrum and 



T h (f) = / dO Ff(n)F?(fl)e 



27rifdQ-s 



is the (non normalized) overlap red uctio n func tion for spin 2 waves. 
The most general expression for (4.11) and (4.14) can be shown to be 



T w (t) = it {AJt) Tr(Z>!) Tr(2? 2 ) + 2B w (t) Tr^ V 2 ) 

+ C w (t) [Tv(V 1 )Tt(SV 2 ) + Tv^TiiSV,)} 

+ AD w (t)Tt(SV 1 V 2 ) + E w {t)Tt{SV 1 )Tt{SV 2 )} 



(4.13) 



(4.14) 



(4.15) 



w 



where r = 2irfd, S = s ® s, and T>k is the tensor of the fcth detector. Following the procedure sketched in Refs. 

§, p9||3C| ] , the coefficients A, B, C, D, and E can be expressed as linear superpositions of Bessel functions (see Appendix 
). The traces appearing in (4.15) carry information about the geometry and the relative orientations of the detectors 
that are correlated. 



B. The noise power spectrum 



d. Resonant mass detectors 

Let us consider a generic multi mode resonant mass antenna, with the modes labelled by an index N, as for example 
a resonant sphere, where N = nl. 

The noise power spectrum is a resonant curve peaked at the proper frequency f N of the modes. It can be char- 
acterized by its value at the peak S n (f N ) and by its half height width, which gives the bandwidth of the resonant 
mode. 

We now generalize the results of Refs. [ plp7| ] concerning S„(f N ) of resonant spheres to the scalar tensor theory. 
Denoting by /3/v the transducer coupling factor (the fraction of the total mode energy available at the transducer 
output), in the case of spin 2 GWs we have 

q ft. * \ 1 4irGhf3 N 47rhp n2 

S n (h;f N ) = 7—2 — ^ = w 3P 1 (4- 16 ) 

l + ai E w Mv z F n 



where the expression (3.15) for E N = S^(n; I = 2) has been used. 

For spin GWs, since we have considered only the first order terms in the expansion in powers of q a of the integrated 
cross section, we write 

= T^% (l-2^)^, (4.17) 



where S^r is now the obvious generalization of (3.25), (3.26) to the case q a ^ 0. Making explicit (4.17) for both the 
monopole and the quadrupole modes, we get 



Finally, let us remark that formulae like ( 1.16 ), ( 4.17 ) hold in general for any kind of detector, because the depen- 



dence of the cross section on the coupling constants G, oP a and q a doesn't change according to the geometrical features 
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of the antenna itself. Actually, if no scalar fields are present, a° a = q a = 0, (4.17) vanishes and (4.16) becomes the 
well known formula (4.5) of Ref. pl[] for the maximum sensitivity to spin 2 waves. 
The bandwidth of the resonant mode is given by 



•/a 



.p-1/2 



(4.20) 



where Q N is the quality factor of the mode, which is of the order of 10 7 and T N is the ratio of the wideband noise to 
the narrowband noise in the JVth resonance mode. 
e. The VIRGO interferometer 

In the frequency region above 2 Hz the noise power spectrum of the VIRGO interferometer can be approximated 
by the following analytical expression [32|p| 



SnU) = Pi 



fo 



P2 



fo 




ith 



/o = 500Hz, Pi =3.46x10" 



'Hz- 



9 x 10 _46 Hz~ 



(4.21) 



3.24 x 10" 46 Hz 



In this parametrization Pi and P-i give the contribution of the pendulum and its internal modes to the thermal noise, 
respectively. P3 controls instead the shot noise contribution. For frequency smaller than 2 Hz we assume that the 
noise power spectrum goes to infinity. 

C. Sensitivity of a pair of resonant spheres 

We consider now the correlation between two resonant spheres. As we are interested in scalar waves, we compute 
the correlation between the monopole modes. Since the monopole tensors are isotropic 



v 



(00) 



1 



(4.22) 



the overlap reduction function depends only on the frequency / and the relative distance d. From the general formula 
(4.15) one finds 



r e (r) = 7T 



^(r) + |iJ f (r) + |c c (r) + Df(r) + \e c (t) 



which explicitly reads (see Appendix |b|) 



r £ (/)=47rj (r). 



(4.23) 



(4.24) 



In Fig. ^ we plot this function for d = 50 km, which is roughly the minimum distance to decorrelate seismic and 
e.m. noises and for d = 400 km, which is the distance between the sites of the resonant bars NAUTILUS, in Frascati, 
and AURIGA in Legnaro. For d = 50 km we observe that the first zero of the function is around 3 kHz, which is a 
frequency higher than the first resonant frequency for both the solid mass and hollow sphere. For d — 400 km, the 
first zero moves back at around 400 Hz. 

Restricting ourselves to metric theories (q a = 0), the SNR is 



SNRc 



'In 



1 + al 8tt 3 



^ T ^c(/»o) r c(/"o) 

2 "° /nO^nte fno) 



(4.25) 



6 With respect to this reference, the value of P2 is slightly changed as can be found in frittp: / /www. virgo.infn.it / senscurve . We 
quote here the most recent value. We thank the referee for pointing this out. 
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FIG. 1. The overlap reduction function F^(/) for two resonant spheres, located at relative distance d = 50 km (solid line) 
and d = 400 km (dashed line). 




where /„o is the resonance frequ ency for the nth monopole mode. To evaluate (4.25) we need the nois e pow er spectrum 
at resonance / = f n o given by ( 4. IS ) where we assume /3„o = 0.1 pl[ . The only free parameter in ( 4.25 ) is a§ which 
can be conveniently expressed in terms of the post Newtonian Eddington parameter fllQfl as 



1 



7Edd 



(4.26) 



7Edd can be measured in light deflection experimen ts |33| ] ■ The minimum detectable scalar spectrum, for an observation 
time T = 1 year, is found imposing SNR^ = 1 in (4.25). To isolate its dependence from 7Edd, in Table || we list the 
reduced spectrum fL 



fijj (1 — 7Edd)/2. We will comment in the Conclusions on the results obtained. 



In Tables ||, H, IV, we consider hollow and solid mass spheres made of CuAl, A15056 and Mo, materials which 
have high density and high velocity of sound |2l],|25| . The geometrical features of such spheres are the outer diam eter 
$ and the ratio £ between the inner and outer radius. We consider the first excited monopole mode, so in ( 4.18 ) we 
put n = 1 and fio = fo- For solid mass spheres we have Hi = 1.14 |B3 | , while for hollow spheres Hi is a function of 
C and some interesting values are listed in the tables themselves pSf T^ 

TABFE II. Minimum detectable scalar spectrum (SNR^ — 1, T — 1 year) for the correlation between the first monopole 
modes of two hollow spheres, with 6 meters outer diameter, made of CuAl (v = 4700 m/s) at d — 50 km. 



M (ton) 


c 


Hi 


fo (Hz) 


A/ (Hz) 


r c (/o) 


y/S n & /o)(Hz-^) 


h 2 C) 

"-0 s '7Edd 


832 


0.25 


0.73727 


770 


24.3 


11.2 


2.21 x 10" 24 


4.5 x 10" 8 


740 


0.50 


0.49429 


609 


19.2 


11.7 


2.86 x 10 -24 


4.0 x 10" 8 


489 


0.75 


0.4307 


498 


15.7 


12.0 


3.77 x 10~ 24 


4.1 x 10" 8 


230 


0.90 


0.42043 


455 


14.4 


12.1 


5.57 x 10~ 24 


7.1 x 10" 8 



TABFE III. The same of Table 0in the case of two 31 ton Mo (v = 5700 m/s) hollow spheres. 



$ (m) 


C 


Hi 


fo (Hz) 


A/o (Hz) 


r e (/o) 


VS n (£; /oXHz- 1 ^) 


h 2 C> 

"0"7Edd 


1.82 


0.25 


0.73727 


3027 


95.7 


0.1 


9.5 x 10" 24 


2.3 x 10" 3 


1.88 


0.50 


0.49429 


2304 


72.9 


3.5 


1.1 x 10 -23 


6.2 x 10 -5 


2.16 


0.75 


0.4307 


1650 


52.2 


7.2 


1.2 x 10 -23 


1.5 x 10 -5 


2.78 


0.90 


0.42043 


1170 


37.0 


9.6 


1.3 x 10" 23 


4.8 x 10" 6 



Similar analysis of correlations can be repeated for the quadrupole vibrational modes of the resonant spheres, which 
can be excited by both spin and spin 2 waves. 
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TABLE IV. The same of Table [n| in the case of two solid spheres made of CuAl («=4700 m/s) and A15056 («=5440 m/s). 





M (ton) 


$ (m) 


fa (Hz) 


A/ (Hz) 


r € (/o) 


y/Sn(£; /o)(Hz~ i/;i ) 


, 2 n 

"0"7Edd 


CuAl 


105 


3 


1672 


52.9 


7.1 


5.0 x 10~ 24 


2.6 x 10" 6 




167 


3.5 


1433 


45.3 


8.3 


4.0 x 10~ 24 


9.2 x 10" 7 




250 


4 


1254 


39.7 


9.2 


3.2 x 10" 24 


4.0 x 10" 7 


A15056 


38 


3 


1935 


61.2 


5.6 


7.2 x 10~ 23 


1 x 10~ 5 




60 


3.5 


1658 


52.4 


7.1 


5.7 x 10" 24 


3.2 x 10" 6 




90 


4 


1451 


45.9 


8.3 


4.7 x 10" 24 


1.3 x 10" 6 



Let us face the calculation of the overlap reduction function. The quadrupole tensors T>^f are traceless, therefore 



the only non vanishing terms in the overlap reduction function, for any (e, e'), are 



r («')( 7 



2B w (t)Tt (v[ e) P^ ) + 4D w (r)Tr (sv{ e) V { f) + E w (t)Ty (sv[ e) ) Tr (s D^) 



(4.27) 



So, the functions to be inserted in the SNR in order to take into account all the cross correlations between the two 
spheres are 



(4.28) 



In FigJ| we compare the plots of Th{f) and The former was computed in Ref. [Q. We choose d = 400 km, 

the distance between the sites of AURIGA and NAUTILUS. 

FIG. 2. Correlation between the quadrupole modes of two spheres, one located at the site of AURIGA (45.35 N, 11.95 E) 
and the other at that of NAUTILUS (41.80 N, 12.67 E), d = 400 km: comparison between F h (f) (solid line), and r f (/) (dashed 
line) . 



Hz 



D. Sensitivity of VIRGO with a resonant sphere 

In this subsection we evaluate the minimum detectable scalar spectrum correlating the monopole mode of a sphere 
and an interferometer, with the noise power spectrum of VIRGO. We label T>i the tensor of the monopole mode of 
the sphere and T>2 the one of the interferometer. Since Tr(Z>i) = 3/2 and T>2 is traceless, the general expression for 
the overlap reduction function is 



3Q(r) + 4D s (r) + J%(t) 



Tr(5P 2 ) . 



(4.29) 
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Explicit evaluation shows the overlap to be 



r e (f)=47rj 2 (r)Tr(SV 2 ) 



(4.30) 



As pointed out in Ref. [^) , this function is a product of a part depending only on the distance d and on the frequency 
and a part depending on the relative position and orientation of the detector frames. For an explicit estimate of the 
minimum detectable scalar spectrum we need t he expression of the interferometer detector tensor with respect to the 
Earth centered reference frame, so that ( 4.30| ) writes in terms of the latitude and longitude of the antennas. The 
interferometer tensor is 



1 



(cos 2\ — cos 2ip) e+ (f ) + (sin 2\ — sin 2tp) e£ (r) 



(4.31) 



where \ an d 4> are the orientations of the two interferometer arms measured counterclockwise from the true North. 
Therefore, the location dependence of this tensor is split into a part depending on the position of the interferometer 
on the Earth surface and a part depending on the orientations of the arms with respect to the true North]]. 

Since the sphere noise power spectrum is narr owba nded with respect to that of an interferometer, we assume the 
latter to be constant, and equal to the value of ( 1.21 ) for / = /„o, within the sphere bandwidth A/„o- This implies 
that the SNR can be written as 



SNR 5 = 



3HZ 



—At T 

2 ln0 



^c(/no)r^(/„o) 



(4.32) 



In Fig. U we plot T^(f) for d — 58 km and d — 270 km, which is the distance between VIRGO and the site of 
NAUTILUS in Frascati. In this figure we plot also the curve for the case in which the sphere is located nearly in 
the Gran Sasso underground laboratory (d = 294 km). For this correlation T^(/) gets approximately its maximum 
value. This result is particularly important in view of the remark that for the future resonant detectors with project 
sensitivity approaching the quantum limit, the cosmic ray interactions in the detector may set a limit to the sensitivity 
in an unshielded environment. 

FIG. 3. The overlap reduction function T 5 of VIRGO (43.63 N, 10.50 E; x = 71.5 deg, ip = 341.5 deg) with a resonant sphere 
located at: Frascati (41.80 N, 12.67 E), d=270 km (solid line); (43.2 N, 10.9 E), d=58 km (dashed line); Gran Sasso laboratory 
(42.4 N, 13.70 E), d=294 km (dotted line). See the text for further explanations. 




For d=270 km the function has its first peak at ~ 591 Hz, with ~ 0.98, and its first zero at ~ 1019 Hz. Notice 
that, as the spacing d increases, the peak frequency moves to lower values. This is peculiar of the correlation between 
a sphere and an interferometer, because for two spheres, two interferometers or two bars an increase of d simply 
implies a shift to lower values of the first zero. For a sphere and an interferometer this effect is due to the function 
j2(r) in (4.30). The relative orientation dependent factor, Tr(6> P2), accounts for the full amplitude of the overlap 
reduction function, being an oscillating function. In fact we can express the overlap reduction function with respect 
to the natural frame of the interferometer defined before H], In this frame, the direction of the unit vector s joining 



7 Notice that the computations of Ref. M were performed in the interferometer frame and not in the Earth centered one. 
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the antennas is determined by the angles (6,4>) and the detector tensor has the simple form (3.2). Using the same 
convention for the angles as in we have 

Tr(5 V 2 ) = ^ sin 2 9 cos 20. (4.33) 

In Table ^ we consider Mo and CuAl hollow spheres at d = 270 km. The CuAl sphere gives the best sensitivities, 
because the resonance frequency is lower and it belongs to the range where the overlap reduction function gets its 
maximum. Sn is the noise power spectrum of the sphere and Sn that of VIRGO, both for the frequency /io = /o- 

TABLE V. Minimum detectable scalar spectrum (SNRj = 1,T = 1 year) obtained by the first monopole vibrational mode 
of one hollow sphere at Frascati with VIRGO (d = 270 km). The CuAl and Mo hollow spheres are those of Table [n] and Table 



[II, therefore only their resonance frequencies and the corresponding sensitivities are written here. 





fo (Hz) 


A/o (Hz) 


r 4 (/o) 




TsFOHz 1 / 2 ) 


h 2 o 

"0"7Edd 


CuAl 


770 


24.3 


0.76 


2.2 x 10~ 24 


4.1 x 10~ 23 


1.2 x 10" 5 




609 


19.2 


0.98 


2.9 x 10~ 24 


3.9 x 10~ 23 


6.6 x 10" 6 




498 


15.7 


0.92 


3.8 x 10~ 24 


3.9 x 10~ 23 


5.6 x 10~ 6 




455 


14.4 


0.85 


5.6 x 10" 24 


4.0 x IO -23 


7.2 x 10" 6 


Mo 


3027 


95.7 


0.19 


9.4 x 10~ 24 


1.1 x 10" 22 


1.8 x 10~ 2 




2304 


72.9 


0.16 


1.2 x 10~ 23 


8.6 x 10~ 23 


1.0 x 10" 2 




1650 


52.2 


0.08 


1.2 x 10~ 23 


6.4 x IO -23 


7.0 x 10" 3 




1170 


37.0 


0.35 


1.3 x 10~ 23 


5.0 x 10~ 23 


5.2 x 10" 4 



V. CONCLUSIONS 



The aim of this paper was to study the detectability of scalar GWs from the cosmic stochastic GW background. 
Before discussing our results, we would briefly like to recall what is the best strategy to perform such a measurement 
in our opinion. As we argued in the introduction, as far as the type of detector to use is concerned, it does not seem 
practical to use only L shaped interferometers. If the impinging monopole mode of a scalar GW moves along the z 
axis and the arms of the interferometer lie in the x, y plane, then the two arms will be stretched by the same quantity 
SL. In this case an interferometer set to work on a dark fringe will not detect any signal. On the contrary, in this very 
configuration a spin 2 GW will give its maximum effect. This is a limiting case though. In general the direction of the 
impinging GW will form a certain angle with z and the perturbation due to a scalar will be tangled with that of the 
spin 2 in an inextricable way. In principle one could reconstruct the directional sensitivity patterns for the two spins 
to separate the two signals. The data needed to do this will require much time to be gathered and for this reason also 
this proposal does not seem practical to us. What would happen if the detector were a sphere? Let us first analyze 
the case of a GW from the viewpoint of a reference frame centered in the origin of the sphere. As it was discussed in 
Ref. the GW interacts with the resonant mass detector through the so called electric tensor Eij — i?iojo- If the 
direction of propagation of the incoming GW is along the z axis, the six components of this tensor can be expressed, 
using a null tetrad, in terms of the so called Newman-Penrose parameters which can also be expressed in the basis 
given by the sf™ 1 ^ defined in Appendix C 1 . In turn this basis can be put in relation with the actual measurements 
performed by the detectors [^2| . In this reference frame there is a direct correspondence between the Newman-Penrose 
parameters and the spin of the incoming GW and the information about the spin of the GW can be easily extracted. 
If the direction of propagation of the GW makes an angle with the z axis, the situation is more complicated as it can 
also be seen from (3.14) which is the cr oss s ection given by a GW with polarization e + , e x and from (3.24), which is 
the analogous with polarization e s . In (3.12), for example, even if the incom ing wave is pure spin 2 all the modes are 
excited and the situation is indistinguishable from that described by ( 3.24 ). The situatio n rese mbles very much that 
for an interferometer. But in the case of the sphere, the pure monopole mode given from (3.20) is only excited by the 
scalar wave, giving a clear signal of the presence of a scalar wave. This is the motivation for our proposal to couple 
an interferometer to a resonant detector of spherical shape. Let us now look at the results we have obtained that, for 
the sake of generality, encompass also the sensitivities for pairs of resonant mass detectors. 
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The behavior of the overlap reduction function for pairs of resonant mass detectors (Fig. |lj-|2j) is quite different from 
the one for a resonant mass detector and an interferometer (Fig. |J). In Fig. we see a constant function which 
abruptly goes to zero for certain values of the frequency. These values of the frequency decrease by increasing the 
distance d at which the two detectors are located. The values of the monopole overlap reduction function and of 
Th(f) are of the same order of magnitude: however, the quadrupole T^(/) is even an order of magnitude smaller at 
low frequencies. 

Quite on the contrary, the overlap reduction function for the pair interferometer-resonant mass detector is different 
from zero only in a certain region which depends from the distance between the antennas and the direction s of the 
sphere with respect to the arms of the interferometers. In Fig. || it is shown that the values of the frequencies at 
which the overlap is maximum are in agreement with the resonant frequencies of the planned detectors || . 

The numerical results concerning the sensitivities are given in Tables §0,0 for pairs of resonant mass detectors 
and in Table [v| for the pair interferometer-resonant mass detector. 

The values given in Tables [il, III show the potential of hollow spheres: going from realistic weights for such detectors 
of the order of the dozens of tons to weights of the order of the hundreds of tons (which are non realistic at the present 
state of the art) there is a gain in sensitivity of two orders of magnitude. Such a gai n could also be achieved going to 



materials with a higher speed of sound propagation [ p5| as can be seen from (4.18) where both M and v appear in 
the denominator but v is squared. 

How do these results compare to those obtained for spin 2 GW? It depends on the value of the scalar amplitude 
and the scalar coupling q a since fit is roughly proportional to the scalar amplitude and coupling while f2/j is roughly 



proportional to 1 + a% (see ( |4 . 1 2[) a nd ( 4.1 3| ) ) . At the moment a$ has been only measured in our solar system and its 



value at la level is « 10~ J |33|. Such a small value is given to the fact that Einsteinian general relativity seems 
to be very well verified. Such a value of ojq would give very little chances to the planned resonant detectors to detect 
a scalar GW background which should be limited by nucleosynthesis to be ||(| 

fegfi*(/)d(ln/) < 1(T 5 . (5.1) 

f>10~ 8 Hz 

We have also to mention Ref. |37|, where an estimate of ag was attempted starting from the same nucleosynthesis 
bound: the result is a weak dependence of «q from distance. Our ignorance on the mechanisms that should give 
Q!q its value (cosmological attractor? supersymmetry breaking?) prevents us from further comments on this point. 
Also concocting strategies with resonant mass detectors made of different materials (to exploit the q a dependence) is 
possible, but probably premature given what we said earlier. We remark however that once operating, resonant mass 
detectors of spherical shape could themselves provide a measure of a§ using binary or collapsing stars as emphasized 
by many authors and more in particular in Ref. f25f| . 

As a final comment we remark that the sensitivity of the pair interferometer-resonant mass detector seems to be a 
couple of orders of magnitude less than that of a pair of resonant mass detectors. The plots we have given show that 
a careful choice of where to locate the detectors can account for up to an order of magnitude in sensitivity. 
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APPENDIX A: CORRELATION FUNCTIONS FOR MANY SCALAR FIELDS 



The SNR for scalar GWs has been computed in Ref. following Ref. p0[ |. Actually, although the generalization 
to single scalar theories (Brans-Dicke) is trivial, the one to multi scalar theory needs the introduction of a very strong 
constraint on the fields themselves: in order to get a formula for the SNR one can state the following lowest order 
condition for the correlation function between the Fourier amplitudes of the scalar fields 



< e*(f, m b (f, n) >= 7o fc Kf - - ft)K(f) 



(Al) 



where K(f) is a real non negative symmetric function. This hypothesis means that the correlation function is the 
same for every pair of scalar fields. This is not the most general situation one can imagine: in fact, because of the 
symmetry a <-> b, one would expect to have n(n + l)/2 distinct correlation functions K a b(f). On the other hand, if 
we consider only one degenerate function K(f), then the framework is exactly the same of Refs. |j|,|30||, and, with the 
same algebra, we get (4.12). 

We reproduce the main steps to express K(f) in terms of the spectrum Straightforward generalization of || 
shows that, for any tensor multi scalar theory, the energy density carried by a GW is 



TOO = Ph + Ps = (l 



32ttG 



< >+8 1 ° ab <ee> 



(A2) 



where the brackets < . . . > stand for integration over a finite region of tridimensional space co nta ining several 
wavelengths. From this formula we recover K{f) as a function of the scalar spectrum Vi s (f). From (Al) one obtains 



<re b >- 32 7r 3 7o afc / dffK(f). 
Jo 



Using the definition of £l s , for non negative /, we get 



/ d Ps 

Pc df 



647T d 



(A3) 



(A4) 



where we used 7g 6 = n. Furthermore fl s = ^ a Qja = nil^a so we can infer from ( A4) 



K(f) = (l + al) 



64^ f *W)' 



(A5) 



where = fl^. 



APPENDIX B: THE OVERLAP REDUCTION FUNCTIONS 



In the following we give the coefficients, introduced in the main text, for the functions T$(f) 



and T h (f) 



B 1 

C 
D 



B 1 

C 
D 



/ T 2 jo(T)-2Tj 1 (r)+j 2 (r) \ 

-r 2 jo{T) +4rji(r) -5j a (r) 

T ji( T ) ~ 5.72(t) 
W jo(r)- 10 Tj 1 (r)+ 35 j 2 (t) J 



/-r 2 jo(r) -2Tj 1 (T)+j 2 (T)T 2 j (T) - 2rji (r) + j 2 (r) ' 
-t 2 Jo (t) - 2rji(r) - 5j 2 (t) 
-t 2 Jo(t)+4tj 1 (t)-5j 2 (t) 
r 2 j (r)- 10^)+ 35 j 2 (r) 



(Bl) 



(B2) 
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APPENDIX C: DETECTOR TENSORS 



1. The sphere mode tensors 



A basis for the pure spherical harmonics is given by the S( lm \ with I = 0,2 |§ 

1 




(CI) 

v \ T l -i / 

The normalization is chosen so that sf™^h l n^ — Yi m . h is the radial unit vector. The vibrational response of a 
spherical detector is usually written in terms of this pure spin basis. Otherwise, following Zhou and Michelson [fl9|| , 
the vibrations of a resonant sphere are more conveniently described as functions of the real quadrupole spherical 

1 /a 

harmonics, in addition to the monopolc spherical harmonic Yoo = (47r) 

Y a = Y 2Q , 

Yic = -j= (Y2-1 - Ka+i) , 

r ls = -^=(y 2 _ 1 + y 2+1 ) , (C2) 
y 2c = -L(y 2 _ 2 + r 2+2 ) , 

^28 = — ^ (3"2-2 - ^2+2) ■ 

A convenient basis for the real spherical harmonics is given by £>(°o) = ^5(00) and 25(e) with £ = 0,1c, Is, 2c, 2s. 
These traceless tensors are defined as 




(C3) 



with 
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From the definition of the real spherical harmonics we have 

X>(°) = _./^5(20) 
V 15 



= — i (5( 2+1 ' + S^) , (C5) 

1.5 



p (2c) /27T / 5 ( 2+2)+5(2 _ 2: 



p(2a) =i y^| ^(2+2) _ ^(2-2) 



2. Explicit expressions of the polarization tensors in the detector frame 



Let us consider now the wave frame (m, n, 0) and the detector frame (x,y,z) defined in (3.5) by introducing the 
rotation matrix 

(cos (f> sin \ 

— cos sin cos cos sin , (C6) 
sin sin — sin cos cos / 

where the angles (0, 0) are defined following the conventions of Forward Jl6|]. The polarization tensors of the GW in 
the antenna frame e B (Cl) are obtained by rotating the ones in the wave frame e B as 

e B ({l) = i? t (fi)e B i?(0); B=x,+,a. (C7) 

The tensors in the wave frame are 

/l 0\ /0 1 0\ /l 0\ 

e+ = I -1 , e x = 1 1, e s = 1 . (C8) 
\0 0/ V Q 0/ Vo 0/ 



From (|C7|), for the spin 2 polarization tensors we get 



, ( 2 (cos 2 - cos 2 0shr 0) (1 + cos 2 0) sin 20 



sin 20 sin < 



e + (f2) = m® m - n ® h = - (1 + cos 2 0) sin 20 2(sin 2 - cos 2 cos 2 0) - sin 20 cos , (C9) 



2 



sin 29 sin — sin 20 cos — 2 sin 



2 , 



C- cos sin 20 cos cos 2(f) cos sin \ 
cos cos 20 cos sin 20 sin sin . (C10) 
cos sin sin sin / 



For the scalar polarization tensor we have 



/2(cos + cos 0sin 0) sin sin 20 



sin 20 sin < 



e s (VL) = rh ® m + n ® n = - sin sin 20 2(sin + cos 2 cos 2 0) sin 20 cos . (Cll) 



2 



sin 20 sin sin 20 cos 2 sin 



3. The Earth centered reference frame 



In Appendix [b] we have listed the coefficients which give the dependence of the overlap reduction function on the 
frequency and on the distance between the antennas. To infer its dependence on the relative orientations of the 
detectors, it is convenient to express the detector tensors in the reference frame of the Earth. We then express the 



23 



detector tensors with respect to a triad of orthogonal unit vectors (x,y,f), where x and y lie on the tangent plane 
and r points along the Earth radius. This triad defines univocally the antenna coordinate system. Given the latitude, 
0, measured in degrees North from the equator and the longitude, $, in degrees East of Greenwich, England, the 
relation of the triad of vectors (x, y, f) with respect to the Cartesian reference frame (X, Y, Z) originated in the center 
of the Earth is 

x = - sin cos $X — sin sin <E>F + cos OZ, 

y = -sin$i" + cos$F, (C12) 
f = cos cos <&X + cos sin $Y + sin OZ . 

A simple example is given by the tensor of an interferometer ]lr| which, in the Earth centered frame, is usually 
written as 

V ij (X,Y) = l(X i X j -Y i Y j ), (C13) 

where X and Y are chosen to point in the detector arms directions. 

Frame dependent expressions of the same kind can also be written for the tensors des crib ing the geometrical features 
of the modes of a resonant sphere. In the Earth centered reference frame, equations ( |C5| ) become 

V S ] = it + ^ - 2ii ^) ' 

D 8 a) = -^(w+M-). ( C14 ) 
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